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Abstract 

Even if quark and lepton masses are not uniquely predicted by the fundamental 
theory, as may be the case in the string theory landscape, nevertheless their pattern 
may reveal features of the underlying theory. We use statistical techniques to show 
that the observed masses appear to be representative of a scale invariant distribution, 
p{m) ^ l/m. If we extend this distribution to include all the Yukawa couplings, we 
show that the resulting CKM matrix elements typically show a hierarchical pattern 
similar to observations. The Jarlskog invariant measuring the amount of CP violation 
is also well reproduced in magnitude. We also apply this framework to neutrinos using 
the seesaw mechanism. The neutrino results are ambiguous, with the observed pat- 
tern being statistically allowed even though the framework does not provide a natural 
explanation for the observed two large mixing angles. Our framework highly favors a 
normal hierarchy of neutrino masses. We also are able to make statistical predictions 
in the neutrino sector when we specialize to situations consistent with the known mass 
differences and two large mixing angles. Within our framework, we show that with 95% 
confidence the presently unmeasured MNS mixing angle sin^is is larger than 0.04 and 
typically of order 0.1. The leptonic Jarlskog invariant is found to be typically of order 
10~^ and the magnitude of the effective Majorana mass niee is typically of order 0.001 
eV. 



1 Introduction 



Of the 28 parameters of the Standard Model (including neutrino masses and mixing) the 
quark and lepton masses appear particularly puzzling. We would expect that the spec- 
trum of masses would exhibit some underlying structure, much as the periodic table and 
the hadron spectrum respectively reveal the dynamics of atoms and elementary particles. 
However, apparently the quark and lepton masses do not show any such pattern. Aside 
from a rough correlation of mass with generation, decades of searching have not revealed 
any significant regularity in the masses. Perhaps there is none. 

In this paper, we explore the possibility that the masses and mixings are not uniquely 
predicted, but are representative of an ensemble of possibilities. The masses and mixing 
then reveal the "weight" or "measure" of the underlying theory. This possibility was first 
proposed by Donoghue in and we develop it further. This is a plausible outcome of the 
string theory landscape 0. In this description, there are very many Standard Model vacua 
with parameters close to those observed, and yet many more with different parameters. 
While we observe only a single ground state, the solution that we observe is representative 
of the ensemble of possible SM vacua that occur in the landscape. In particular, the Yukawa 
couplings would not be unique, but would be representative of the couplings found in the 
ensemble of solutions. Since we have many manifestations of the Yukawa couplings (9 
quark and lepton masses, two neutrino mass differences as well as the CKM and MNS weak 
mixing elements) we may apply statistical methods even though we live in a single member 
of the possible landscape vacua. Our primary assumption is that the observed masses are 
representative of this ensemble. By studying the phenomenology of the masses, we then 
learn about the nature of the underlying theory. 

By employing statistical tests we will show that the observed quark masses appear to 
be distributed in a scale-invariant fashion. This was already suggested in but we provide 
better statistical measures. This result can be readily seen in Fig. ^ A scale invariant set 
of masses is one where the values appear as a uniform random distribution when plotted 
on a log scale. The masses of the quarks and charged leptons are shown on such a scale in 
Fig. Hand appear visually to be consistent with being uniformly random. In Section 3, we 
use statistical tests to confirm this and to quantify how far the weight could deviate from 
a scale invariant form. In Section 4 we apply this distribution to the full Yukawa matrices 
of quarks and investigate the weak mixing angles. We will see that one naturally generates 
a hierarchy in the CKM elements similar to the one found in nature, aside from a discrete 
ambiguity concerning the generation structure. Moreover, the magnitude of the observed 
CP violation is readily reproduced. 

The neutrino masses appear quite different from other lepton masses, with a possible 
explanation being the seesaw mechanism and the presence of heavy right handed Majorana 
mass. We explore the neutrino masses in a similar fashion, using the same weight for the 
Dirac masses and allowing a different weight for the Majorana masses. We describe the 
neutrino phenomenology in Section 5. 

It is possible that these results could be useful in testing or exploring the string land- 
scape. For example, it has been suggested jSj that a scale invariant weight could emerge 
from the intersecting brane worlds construction 1| of the Standard Model. If some branches 
of the landscape yield weights compatible with observations while others do not, this could 
be used to refine our location in the space of string theory solutions. We are clearly far from 
an understanding of the landscape, but the masses and mixings of the leptons can play a 
useful role in landscape phenomenology. 
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Figure 1: Quark and lepton masses, defined at tiie energy fj, = Mw, on a logarithmic scale. A scale 
invariant weight corresponds to a uniform distribution on this scale. 

2 The properties of the weight 

The most fundamental assumption of our approach is the existence of a large ensemble of 
nearly equivalent vacuum states with different values of the quark masses, of which our 
vacua is assumed to be representative. The string landscape could be such a theory. It 
has been estimated j2I that there are perhaps 10^'''' Standard Model vacua with parameters 
agreeing with observation within the experimental error bars. There are then many more 
that look like the Standard Model but have other values of the quark and lepton masses. 
This larger set of SM vacua then constitutes our ensemble. 

Within this ensemble we need to make other assumptions. For example, we are assuming 
that the quark masses are not correlated with each other. If there are really of order 10^'^*' 
vacua this seems reasonable. We could hold any one of the quark masses fixed and still find 
a whole range of values for the other masses. We also will make the assumption that quark 
and charged lepton masses (but not neutrino masses) play by the same rules, although 
we will explore some slight differences that arise due to the differences between quarks and 
leptons. Finally, we need to make assumptions about the whole matrix of Yukawa couplings, 
as will be discussed below. 

The weight p{m) is a probability distribution function for the masses (or later for the 
Yukawa couplings). It is defined by considering the fraction of masses / found at a value 
m within dm as 

/(m) = p{m)dm. (1) 
The normalization of the weight is by definition then 

1 = y p{m)dm. (2) 

For the simple weights p[m) that we consider, this normalization condition demands upper 
limits and, in some cases, lower limits for the range of allowed masses. 

Because the masses depend on the energy scale, the weight will also be depend on this 
scale. The renormalization group allows us to determine the scale dependence, as was shown 
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in When the energy scale is changed from //i to fi2 a mass mi will change to the value 
7712, defining a functional relationship m2 = m2(mi) or the reverse mi = 1711(1712). In order 
to preserve the definition of the weight, it must transform as 



with 



dmi 



J 



As an example, under QCD the masses transform as 



with 



m{ii2) = m{iii) 



-I dr, 
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In this case the renormalization group transformation rule yields 
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(3) 
(4) 

(5) 
(6) 

(7) 



For large values of the masses we also need to consider the effect of the Yukawa inter- 
action on its own running. Recall that in the Standard Model the masses rrij are related 
to the Yukawa couplings hi via mi = hiA= where v is the Higgs vev. The renormalization 

v2 

group equations for the combined QCD gauge and Yukawa running for Nf = 6 are 
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dh^ 
dt 
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where t = log(//f //i^)^ using /ii as the initial scale. Note that the QCD interaction tends to 
make the Yukawa coupling larger as one scales down in energy (t positive) , while the Yukawa 
self interaction tends to make the Yukawa coupling smaller. One of the consequences of this 
is the well known quasi-fixed-point 5 at m^, = 220 GeV - all large Yukawa couplings will 
run to a value close to this fixed point when scaled down from large energy to the W scale. 
The solution to the Yukawa RGE has the form j2j 
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From this we can extract the Jacobian for the weight 

52 



Jim) 



62 - 2a'- 



3/2 ■ 



(9) 



(10) 



3 



For small masses, this is equivalent to the gauge rescaling given above. For leptons, we 
neglect the QED effects and only consider the Yukawa interaction. In this case, the above 
formulas continue to hold, but with the identification 

m = 1 

a{t) = (11) 

Let us look at a particular weight which plays an important role in the rest of this paper. 
This is a weight that has the form p(m) ~ 1/m, which we will call the scale invariant weight. 
To normalize this form we need to limit the range of masses with an upper and a lower 
cutoff, so that the complete form is 

1 1 

p (m) = @ (m — m_) B (m-|_ — m) . (12) 

log — — m 

This distribution is scale invariant in two senses. First, a probability distribution that goes 
as dm/m is clearly invariant under a rescaling of all masses m — > Am. The endpoints are 
not invariant but they rescale in the obvious fashion. In addition, this weight is invariant 
under renormalization group transformation of the gauge interactions changing from one 
scale /ii to another scale This can be seen simply from Eq. 0. 

However, the Yukawa interactions modify the shape of this weight. In particular, if the 
weight takes the above form at a scale fii, then at a lower scale /i it will have the form 



where 



1 



2ci Tn — T 



Q (m — m_) Q (m-|- — m) (13) 



g m 



m± = z T7TJ2 (14) 



are the rescaled endpoints. 

We can illustrate the renormalizaton group effects by considering the transformation 
from a Grand Unified scale of /i = 10^^ GeV down to the W scale. At the GUT scale we 
postulate for the Yukawa couplings the weight 

P>^ (h) = lQ{h-h^)@{h+-h) (15) 

log^ h 

with /i+ = 1 and h- = 1.2 x 10^^. The values of the endpoints have been chosen to 
allow the mass ranges to extend from below the electron mass to above the top quark 
mass. When transformed down to the W scale, the maximum value of the quark mass 
distribution is m+ = 197 GeV and the minimum value is m_ = 0.53 MeV. For the leptons, 
the corresponding values are m+ = 103 GeV and m_ = 0.20 MeV. Note that the rescaled 
ranges of the distributions are different for quarks and leptons. At low values of the mass the 
weight retains the 1/m form, but the larger values the shape of the distribution is modified. 
This is shown in Fig. |2l(a) for the quarks and Fig. 121 (b) for the leptons. The modification to 
the shape is gentle enough that we will ignore it in most of the phenomenological applications 
below. However, in the next section we will briefly explore the effect of rescaling the quark 
and lepton masses to the GUT scale and fitting the weight at that scale. 
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Figure 2: (a) The solid curve corresponds to a pure scale invariant weight with an maximum Yukawa 
coupling of 1 and a minimum of 1.2 x 10^^. The dashed curve corresponds to using such a weight 
defined at the GUT scale and evolving it down to the W scale using the renormalization group 
equations for quarks. Note the slight distortion of the shape as well as the shift in the endpoints of 
the distribution, (b) Same as before, except that the renormalization group for leptons was applied. 

3 Quark and charged lepton masses 

In this section we will treat the quark and charged lepton masses and will determine the 
weight that best describes their distribution. We describe the masses at the scale of M]y. 
Running the quark masses to M^y the dominant contributions come from QCD running 
and we can neglect electroweak running for our purposes. The quark masses at the scale 
Mw ^ are given in Tabled For leptons we neglect the small effects on the masses from 
electroweak running so that the lepton masses we use at Mw are simply the ones quoted 
by the PDG as given in Table El 



ruu 


rud 


rris 




nib 


nit 


2.2 MeV 


4.4 MeV 


80 MeV 


0.81 GeV 


3.1 GeV 


170 GeV 



Table 1: Quark masses at the scale Mw- 



me 




m-r- 


0.511 MeV 


0.106 GeV 


1.78 GeV 



Table 2: Lepton masses at the scale Mw- 



From these values, we can deduce several properties. The masses span a range of more 
than five orders of magnitude, from nie to mt- Most obviously, the weight is not flat, but it 
must be peaked towards low masses. It is also not likely to be a Gaussian or an exponential, 
because such distributions would have an exponentially small probability of accommodating 
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the top quark. The weight must have a peak at low mass, yet extend out to high mass in 
order to explain the top quark. 

We will explore the class of power law weights, and will see that they are successful in 
describing the distribution in masses. The power law weights have the form 

p(m)~^. (16) 

Obviously they are not normalizable if the range includes all values of m, so we need to limit 
the range of the values of masses m. For 5 < 1 we need an upper bound of the range of m 
so that the weight is normalizable. A natural choice for that is the quasi fixed point of the 
SM, m* = 220 GeV, which we will use as an upper bound for masses for all values of 6. For 
5 > 1 we also need a lower bound for the range of m. In order to explain the existence of 
the electron but no lighter fermions, a sensible value for the lower bound should be smaller 
than m-e but not by many orders of magnitude. We choose mio^ = 0.4 • where me is 

the electron mass. Varying this value does not significantly affect the conclusions presented 
in this section. Using these bounds and normalizing the weight, we obtain the correct 
expressions for power law weights: 

A) for 5 = 1: 

P [m) = — G (m - miou,) S (m* - m) (17) 

B) for (5 7^ 1 with a lower bound miow- 

P("^) = —r^ ^—r-i^&im-miou,)<d{m^-m) (18) 



C) for 5 < 1 without a lower bound miow'- 

The distribution of a sample of random numbers with a power law weight p (m) ~ 
looks like a uniformly distributed sample on the scale m}~^ (for 5 ^ \). For the scale 
invariant weight p{m) ~ a sample looks uniformly distributed on the scale logm. Fig. 
H shows the physical masses of the quarks and leptons on a log scale and it is visible that 
these are consistent with a uniform distribution. For contrast, in Fig. I^Jwe plot the quark 
and lepton masses on a ^/m scale which corresponds to a weight p{m) ~ It is evident 
that this sample is quite unlikely to result from a uniform random distribution. 

We can make this mathematically precise by using the maximum likelihood method to 
select the optimum power 5. The likelihood function is defined as the product of the weights 
of the masses: 

L{5) = p{mi,6) for quarks only (20) 

i=u,d,s,c,b,t 

^i^) — pi^ijS) for quarks and leptons (21) 

i=e,fi,T,u,d,s,c,b,t 

The optimum power delta is where L(6) has a maximum Lmax, which coincides with the 
maximum of the log likelihood function log L{6). The one sigma range is limited by the 
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Figure 3: Quark and lepton masses, defined at the energy fj, = AIw, plotted versus \/fri. A weight 
p ~ Xj ^pm corresponds to a uniform distribution on this scale. 

values of delta where logL(5) = logLmax — \ and the two sigma range is limited by the 
values of delta where log L{S) = log Lmax — 2. 

We do this first in the range of masses where there exists a lower cutoff miow which we 
take at 0.4 of the mass of the electron. In this case all powers of 6 are allowed in principle. 
We perform the analysis for only the quark masses and for the quark and lepton masses 
combined. The log likelihood function log L {5) for the quarks and leptons combined is 
plotted in Fig. [l] In this case the favored value of 5 is almost exactly unity, with a range 

<5 = 1.02 ± 0.08. (22) 

If we only use the information in the quark masses we find nearly the identical result, 
5 = 0.99 lb 0.10. The error bars in these results include the effects of the limited statistics 
and get smaller when we include the lepton masses. 

Alternatively we could consider power law weights that extend all the way down to 
m = 0. In this case the power 6 must be strictly less than unity. We obtain 6 = O.SGIq'qj 
for quarks and leptons and 6 = O-SS^'^qq^ for quarks only. All results from our likelihood 
fits are graphically summarized in Fig. El 

Another way to probe the character of the weight is by a moment analysis. A random 
variable x which is uniformly distributed between and 1 has moments 

(x") = (23) 
^ ' n+1 ^ ^ 

Of course for a small number of values drawn from a random distribution the moments 
will not be precisely these values, but will exhibit a scatter around these values due to the 
limited statistics. We can understand this effect by simulations. As an example we consider 
the second moment {x'^) of six random variables drawn from a uniform distribution from 
zero to one. We simulate six random numbers, evaluate {x"^) and repeat that many times. 
That results in a distribution of (x^) which has the shape shown in Fig. |2l From this figure 
one can not only verify the expected average, but also understand the la and 2a ranges. 
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Figure 5: Results of the likelihood method. The blaek ranges are 2(t excluded, the grey ones are 
la excluded and the lines inside the white regions are at the optimum powers S. The upper results 
are the ones taking only the quark masses into account, the ones in the middle are from quarks and 
leptons combined and the results at the bottom are the ones for quarks and leptons with anthropic 
constraints. 

For example, at two standard deviations we find the allowed range of the second moment 
of six randomly selected numbers is 

0.112 < (x^) < 0.585. (24) 

We apply these tests to the weights, asking if the masses at Mw correspond to a uni- 
formly distributed random variable in m^^^ or logm. We use the first three moments 
n = 1,2,3 and require that all of these moments fall within the la or 2a ranges. Besides 
the moments (a:"), we also take into account ( (1 — a;)" ) and ( (5 — x)^ ) fov n = 1, 2, 3 
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Figure 6: Distribution of (a:^) for six random variables. The dashed hnes mark the lu ranges and 
the solid lines mark the 2a ranges. 

and always take the most stringent bounds of all of them. The results from this moment 
analysis are plotted in Fig. [7J 

Comparing with the la and 2a ranges from the likelihood method, we find that the 
allowed ranges are somewhat smaller from the moments analysis but otherwise are in good 
agreement. The fact that we always take the most stringent bounds in the moment anal- 
ysis probably lets the allowed ranges shrink in comparison to the ones from the likelihood 
method. 

One subtlety we did not discuss so far is the distortion of our analysis by possible an- 
thropic constraints. This is a difficult subject and we cannot fully address it here. However, 
anthropic issues have the potential to undermine our fundamental assumption that the 
masses which we see are representative of the Standard Model ensemble. This is because 
the formation of complex elements which are necessary for life can only occur for certain 
configurations of the charged fermion masses. In particular some masses need to be light 
jni 13 IB] on the QCD scale - otherwise the nucleons would decay rather than be bound in 
the nucleus. Thus the masses that we see may represent an anthropic bias in favor of light 
masses. 

We can partially address anthropic issues by the following strategy. We will eliminate 
from consideration the up and down quarks and the electron, and will revisit our analysis 
only considering the second and third generations of fermions. This is because the first 
generation masses may be biased by anthropic considerations while there are no known 
anthropic constraints on the heavier fermions. 

Specifically, we proceed as follows. We drop the light masses m„, m^^ and me from our 
analysis. We also exclude the possibility that more than these three fermions are very light. 
This is not a firm requirement as there could be different pattern of elements if more quarks 
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Figure 7: Results of the moment analysis. The black ranges are 2a excluded and the grey ones are 
Icr excluded. The upper results are the ones taking only the quark masses into account, the ones 
in the middle are from quarks and leptons combined and the results at the bottom are the ones for 
quarks and leptons with anthropic constraints. 



or leptons were light. However, this assumption that only the three members of the first 
generation play a role in the elements does correspond to the anthropic setting that we find 
ourselves in, so it makes sense to study only this situation. We apply this constraint by 
not allowing any other masses to occur below 10 MeV. This boundary is chosen because 
of the nuclear binding energy of 10 MeV/nucleon. Fermions with masses greater than this 
do not play a role in nuclear binding because they decay weakly rather than being stable 
in the nucleus. For distributions with a low mass cutoff, we accomplish this by taking 
mioui = 10 MeV. 

The results from our likelihood and moment analysis with anthropic considerations are 
included in Fig. ElandEl It is interesting that this modification does not significantly change 
our results. For weights with a lower bound m^o^ on the masses all results are consistent 
with a scale invariant weight 5 = 1. For the weights extending all the way down to m = 0, 
a power 5 in the range 6 = 0.7 — 0.9 seems preferred. One can understand why there is 
little change in our results by looking at the visual representation of the distribution in Fig. 
n Our procedure for simulating anthropic constraints consists of dropping the three lowest 
masses and applying a cutoff at 10 MeV, which on this figure is just above the mass of the 
down quark. One can see that the remaining 4 quark and 2 lepton masses are consistent 
with uniform distribution by themselves. Thus we have evidence for the scale invariant form 
even considering only the heavier masses for which there is no known anthropic constraints. 

Finally, let us now consider the effect of scaling the quark and lepton masses up to 
the GUT scale (taken to be 10^^ GeV). We could equivalently scale down a distribution 
from the GUT scale and apply separate weights for quarks and leptons, as illustrated in 
Section 2. However, it is simpler to use the renormalization group to scale the masses up 
and study a common distribution at the GUT scale. This explores the issue of whether the 
quarks and leptons should be combined using the same distribution at the W scale or at 
the GUT scale. We will see that there is very little difference between these alternatives. 
The renormalization group equations were given in the previous section. We use only the 
Standard Model interactions in the running of the masses up to the GUT scale. The 
resulting Yukawa couplings defined at 10^^ GeV are given in Table |21 

Again we can visually explore the the comparison of these results with the scale invariant 
weight by displaying them on a log scale, as in Fig. |H1 While there is a small shift in the 
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hu 


hd 


hs 


he 


hb 


ht 


0.46 X 10"^ 


0.92 X 10"^ 


1.7 X 10-4 


1.7 X 10"^ 


6.5 X 10"^ 


0.63 



he 


h^ 




2.9 X 10-*^ 


6.1 X lO""* 


0.01 



Table 3: Quark and lepton Yukawa couplings at the scale Mqut = 10^® GcV. 




Figure 8: Quark and lepton Yukawa couplings, defined at the GUT scale ji — 10^^ GeV, on a 
logarithmic scale. This differs from Fig. Hby the (different) renormalization of the quark and lepton 
couplings as one transforms them to the GUT scale. 



relative placement of the quarks and leptons, due to their different running up to the GUT 
scale, this make no qualitative difference in the apparent randomness of the distribution. 
This is also manifest in a maximum likelihood fit. We consider the weights with both upper 
and lower cutoffs - in particular we use /i+ = 1 and /i_ = 1.2 x 10~^. As noted previously, 
these correspond to a maximum quark mass of 197 GeV and a minimum lepton mass of 
0.20 MeV. A maximum likelihood analysis of the weight for the combined quark and lepton 
Yukawa coupling as a function of 5 again favors values of 5 near one with a fit value of 
5 = 1.06 lb 0.08. The 2a allowed range of 5 is 0.89 — 1.24. These ranges are very similar 
to the ones that we obtained when we simply combined quark and lepton masses at the W 
scale. We conclude from this that effects of different treatments of quarks and leptons are 
small. Also we have learned that it does not make any practical difference if we determine 
the weight at the W scale or at the GUT scale - the distribution is close to scale invariant 
at either energy. 



4 Quark mixing 

In this section we extend our discussion to the full Yukawa matrices for up- and down-type 
quarks, generating not only the quark masses but also the CKM mixing elements. Again, the 
assumption is that each of the Yukawa elements can be treated as an independent variable 
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to be generated independently from a probability distribution function. Of course, even in 
the landscape, this assumption could be incorrect if there are symmetries or dynamics that 
constrain the independent variables in the Yukawa matrices. However, we will see that this 
ansatz does statistically account for the hierarchy that we see in the CKM elements. The 
reason this occurs is simple - because the weight is peaked at low mass, many of the Yukawa 
couplings are small and therefore generate small mixing angles. This dominance of small 
Yukawa elements may also provide an approximate explanation for the idea of "textures" 
in the quark mass matrices [TP. 

The generation of the quark mass matrices is a well-known procedure in the Standard 
Model The mass matrix for the up-type quarks, for example, follows from the Yukawa 
coupling to the Higgs through the relation 



< = (25) 




where hij are the Yukawa couplings. This consists of 9 independent complex elements. The 
diagonalization of this matrix occurs through separate transformations on the left handed 
and right handed quarks 

=1 O" me I = y^")^M^"V^"\ (26) 



A similar process occurs for the down-type quarks 

/ rud \ 
mW= = y^'^^^M^'^Vi'). (27) 

V nib / 

The CKM matrix is then formed by the product of the two left handed rotation matrices 

VcKM = vt^^vt^. (28) 

In general the elements of the mass matrix M^ij are complex- valued. In our model the 
magnitude of the elements of Moij are being distributed from miow to according to 
the weight We take the phases of the elements of Moij to be uniformly distributed 

between and 27r. The appearance of complex elements for Mq^j is required for the existence 
of CP violation and of course is not forbidden by any of the symmetries of the Standard 
Model. However, we have also simulated the mixing matrices with real elements and we note 
that this does not change any of our qualitative results, aside of course from the measure 
of CP violation. 

A first question to be addressed is whether an assumed weight for the Yukawa couplings 
will lead to the same weight being applicable for the masses that emerge as output from the 
diagonalization process. We cannot address this question analytically, so we approach it 
via simulations. The result that we find is that there is a change in the weight between the 
initial Yukawa couplings and the final masses and that this change is roughly categorized 
by a shift in the power 6 by about —0.16. 

In order to demonstrate this, we have generated a large sample of mass matrices of 
Yukawa couplings with a scale invariant weight p ~ l//i, diagonalized the matrices and bin 
and plot the resulting mass values on a logarithmic scale in Fig. |21 (a). If the resulting 
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Figure 9: Elements of mass matrices liave been distributed according to the scale invariant weight. 
In (a) we bin and plot the the eigenvalues of the mass matrices on a logarithmic scale, whereas in 
(b) we do the same on an mP'^^ scale. The vertical dashed lines indicate the values of miow and m*. 

distribution was also scale invariant, this plot would be flat. However, the result is not flat. 
However, when plotted versus mP'^^ in Fig. 111(b), we see that the resulting distribution is 
almost flat, indicating that the result is compatible with a weight of /? ~ Likewise 
if we simulate a weight for the Yukawa interactions of p ~ l//i^'^^ and bin and plot the 
resulting masses on a logarithmic scale in Fig. EH we see an almost flat distribution which 
would correspond to a scale invariant weight for the output masses. 

Another outcome we see is that after the diagonalization of the Yukawa matrices the 
limits miow and are washed out for the mass eigenvalues. There are masses smaller than 
miow and larger than m*, but as one can see in Fig. l9l and [1(11 masses which are orders of 
magnitude smaller than m/o^ or larger than m,, are highly suppressed. 

A weight for the Yukawa matrix elements that is capable of explaining the observed 
quark masses must therefore have a higher power 5 than what we found in our analysis in 
Section 3. As we saw, 6 = 1.16 for the Yukawa matrix elements roughly corresponds to a 
scale invariant distribution for the masses. The previous limits miow and can also be 
used for the Yukawa matrix elements. 

Now we turn to the details of the simulation of the CKM matrices. First, we generate a 
large sample of mass matrices for both up-type quarks and down-type quarks. Each element 
of a particular mass matrix is taken to be independent of each other. By convention the 
largest masses in the two mass matrices are the top and the bottom, so we use the freedom 
to relabel the fields in order to place the largest element in the (3,3) position of the mass 
matrices. 

As discussed above, the diagonalization of the mass matrices is being done by biunitary 
transformations. In general, in this procedure the physical mass eigenvalues of quarks have 
no constraint to fall within any particular range. However for this mixing analysis we wish 
to consider only output masses that are somewhat similar to those observed in nature. 
We impose this by considering only those mass matrices for which mass eigenvalues for 
up and down quarks {i.e the lightest eigenvalues of each matrix) are less than 10 MeV, 
mass eigenvalues for charm and strange quarks are within 50 MeV and 2 GeV and mass 
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Figure 10: Distribution of the eigenvalues of the mass matrices, on a logarithmic scale, when the 
elements of the mass matrices have been distributed according to weight p ~ 

eigenvalues for top and bottom quarks are greater than 3 GeV. 

We construct CKM elements from two left handed rotation matrices, one from up-type 
quarks and the other from the down- type quarks. In the results, there are a few discrete 
classes for the pattern of mixing angles. Since we have not imposed the generation structure, 
there is a set of cases where the up and the strange are dominantly coupled to each other 
and Vud is small. Likewise there is a small number of cases where the charm and bottom 
quarks are dominantly coupled to each other. We discard these cases and only consider 
those results which exhibit the proper generation structure such that the largest couplings 
are to members of the same generation (as defined by the relative masses). 

We use a weight for the Yukawa couplings of which generates the almost scale 

invariant weight for the masses. With these constraints, the resulting distribution of the 
magnitude of the Vus elements of CKM matrices is shown in Fig. 1111 One can see that the 
distribution is peaked at lower values. When one looks more closely at the distribution of 
Vus at very low values, one sees that the peak is located at about 0.002.^ Another feature 
of the distribution is that it vanishes for values larger than l/\/2. The reason for that is the 
generation structure we impose where we require the diagonal elements to be the largest 
elements in their row and column in the CMK matrix. 

^When doing the same simulations without miow this peak does not occur but the distribution for Vus 
reaches a maximum at zero. One can understand that from dimensional analysis: in the simulations with 
"T-ioiD there are two dimensionful parameters involved, m« and m;om where miom/'m* <C 1. However, in the 
case without miow there is only one dimensionful parameter. Therefore a peak for a dimensionless quantity 
such as Vus is expected to be found either at zero or at a value of order one in a case where only one 
dimensionful parameter exists. 
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Figure 11: Distribution of the simulated values of Vus- 

The other CKM elements we consider are Vub and Vcb- For both we observe shapes of 
their distributions similar to the one for Vus as seen in Fig. 1111 but with a higher preference 
for small values. 

For the distributions of the CKM elements one cannot define la and 2a ranges as for 
Gaussian distributions since the peaks are very close to zero but the tails of the distributions 
are extremely long. Therefore we instead consider the medians of the distributions and the 
values where 68% and 95% of the elements lie below them. We call the ranges up to these 
values the la and 2a ranges respectively. For 5 = 1.16 the observed values of Vub and Vcb 
fall within the la ranges and the physical value for Vus is within the 2a range. 

The medians of the CKM elements follow the similar hierarchy that is observed in 
nature. In Fig. we plot the medians of Vub,Vcb,Vus as a function of 6. In general, 
there is a hierarchy in the magnitudes of the elements. This is understandable from the 
diagonalization procedure. Because of the 1/h^ distribution, most Yukawa couplings are 
small. The mixing angles needed to diagonalize the mass matrices are proportional to the 
off-diagonal elements divided by the difference in masses. The mixing involving heavy top 
and bottom quarks will be smaller than those that involve the mixing of down and strange 
quarks. Thus we see that our model can accommodate the hierarchical structure of CKM 
elements. 

After absorbing several complex phases into the definition of the quark fields, the CKM 
matrix can be expressed by four parameters, of which three are mixing angles and one is a 
complex phase. This complex phase indicates the existence of CP violation in the theory. 
Perhaps the best way to describe CP violation is through the "Jarlskog invariant", J \12\ 
which is invariant under rephasing of the quark fields. All CP violation is proportional to 
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Figure 12: Medians of the CKM elements as a function of S 
the product of CKM elements 

J = lm{VudV,sV:,V:,). (29) 

In nature, this is observed to have the value J = (2.88 ± 0.33) x 10"^ In Fig. UHl 

we plot the distribution of J and indicate the la range. We note that observed value of 
J in nature is well within the la range. We conclude that, within the hypothesis of the 
weight, the observed magnitude of CP violation is quite natural for the same probability 
distributions that describe the quark masses. 

5 Neutrino masses and mixing 

The neutrino masses are different enough from those of the other fermions that it is clear 
that they should not be governed by the same weight. Moreover, there are strong theoretical 
reasons why neutrino masses should be treated differently. In this section we explore the 
description of the neutrino masses and mixing in the case where they are described by the 
Type I seesaw mechanism llSl. Specifically we will assume that the Dirac component of 
neutrino masses is governed by the same weight as the charged leptons, and will add the 
possibility of a different weight for the Majorana component of the mass. Unfortunately our 
exploration will be somewhat inconclusive. The observed pattern of masses and mixings is 
not highly typical of the theoretical distributions, but it is also not statistically excluded 
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Figure 13: The Jarlskog invariant for quarks describing the magnitude of CP violation 



at the 2a level. We note that other attempts at simulating neutrino properties with a 
statistical distribution can be found in (18j and ^Hl- We compare our analysis to these at 
the end of this section. 

Our present knowledge of neutrino masses comes from various neutrino oscillation ex- 
periments and from the measurements of the anisotropics in the CMB by WMAP and 
large scale structure formation. From neutrino oscillation experiments, we know two mass 
differences and WMAP gives an upper bound on the sum of the neutrino masses |14| I15j: 

5.4 X 10"^ eV^ < Am?2 < 9.5 x 10"^ eV^ 

1.2 X 10"^ eV^ < |Am?o| < 4.8 x 10"^ eV^ 



^mfc<0.7eV (30) 

k 

First we consider the neutrino masses by themselves without mixing. In the Type I 
seesaw mechanism, the light neutrino masses rrii are of the form 

2 

mi) ,■ 

where mD,i are the Dirac masses connecting left and right handed fields and Mm is the 
Majorana mass term of the right handed neutrinos. 

We assume that the Dirac masses come from the same weight as the quark and charged 
lepton masses p ~ l/m^ and that the Majorana mass Mm is a common scale. Since in 
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the Type I seesaw, the Dirac masses are proportional to the Higgs vev, it is a reasonable 
assumption that the weight is the same as that for the quarks and charged leptons. Note 
that a scale invariant weight for the Dirac masses predicts that the neutrinos will also be 
distributed with a scale invariant weight. If the Dirac masses are distributed with a weight 
~ and the Majorana mass is taken as a constant, then the neutrino masses are 

distributed with a weight 



and the use of 5 = 1 implies 5y = 1 also. 

An initial result is that this scheme greatly favors a "normal" hierarchy in the neutrino 
masses with mi ^ m2 ^ m-3. In this case the solar mass difference is representative of the 
mass scale of the second lightest neutrino with mass m2- The third neutrino is then heavier, 
with a mass ms of order the atmospheric mass difference. The scale invariant weight favors 
this because it favors small masses. It would be unlikely that there are two large masses 
that just happened to be so close together to yield the small solar mass difference. In 
the simulations that we describe below for the case involving mixing, the probability of 
obtaining the mass differences with an inverted hierarchy compared to a normal hierarchy 
was about 0.1%. A quasi degenerate neutrino mass spectrum is even much more unlikely. 
From now on we will assume that a "normal" hierarchy exists. 

One might hope that the observed masses, combined with the assumption that the 
Dirac masses contribute to the seesaw mechanism with the weight determined for the other 
fermions, would lead to a prediction for the scale of the Majorana mass. In practice this 
does not occur. We have performed a likelihood fit for the scale of the Majorana mass 
Mm and we find that the 2a allowed range for Mm spans from 10^ GeV to 10^^ GeV. The 
details of this fit are given in the Appendix. 

One way to see why the limits on the Majorana mass scale are that weak is to plot the 
masses on a log scale. Of course we do not know the absolute masses precisely since we have 
only the mass differences. However, under the assumption of a normal hierarchy we know 
the largest mass quite well and, if the lightest mass is not anomalously close to the second 
lightest mass, we have a reasonable estimate of the second lightest mass. However, there is 
no constraint at all on the mass of the lightest neutrino. Thus the situation is close to that 
pictured in Fig. 1141 where the rough positions of the two heaviest neutrinos are pictured 
and the lightest is allowed at almost any location smaller than these on the logarithmic 
scale. 

In the figure also are shown the ranges allowed for the neutrino masses in the situation 
where the Dirac mass is described by the scale invariant weight and the Majorana mass scale 
Mm is a constant. The ranges are shown for various values of the Majorana mass scale 
Mm- The scale invariant distribution would correspond to a uniform random distribution 
on this logarithmic scale. It is clear visibly that if we allow the lightest neutrino mass to be 
in a broad possible range, that a wide range of Majorana masses is consistent with the scale 
invariant mass distribution. For high Majorana scales Mm there is however more "phase 
space" for the lightest mass which leads to the slight preference for high Majorana mass 
scales Mm as seen in our likelihood fit. Independent of the scale Mm one can see that the 
two heaviest neutrino masses are always very close together on a logarithmical scale. We 
have to consider that a statistical accident in our scheme. 

The analysis of neutrino mixing angles also is not very definitive. Neutrino mixing is 
somewhat different than quark mixing. Let us review the experimental data and the mixing 
formalism within the context of the seesaw mechanism |16) . 




(32) 
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Figure 14: The neutrino masses on a logarithmic scale. The mass differences were trans- 
lated into rough mass estimates as described in the text. Note that the lightest mass is 
unconstrained at small values of the mass. If the distribution is to be described by a scale 
invariant weight, these masses should be uniformly distributed on this log scale. The range 
of the distribution is also shown for various values of the Majorana mass. Given the limited 
statistics, it is hard to determine if a scale invariant distribution is appropriate or to bound 
the Majorana mass scale. 

The lepton equivalent of the CKM matrix is called the MNS matrix. We can parame- 
terize the MNS matrix in the form |10, 



C12C13 



Vmns'- 



-SI2C23 - Cl2S23Sl3e 
S12S23 - Cl2C23Sl3e*'' 



iS 



S12C13 
C12C23 - Sl2S23Sl3e* 
-C12S23 - Sl2C23Sl3e 



-iS 

S23C13 
C23C13 



V 



diagfe*"i/2^e^°2/2,l 



(33) 

where ctj = cos 9ij and Sij = sin 0ij. The three Oij^s are mixing angles and 6, ai, are CP 
violating phases. The phases are only observable if neutrinos are Majorana particles. 
The present 3cj allowed ranges for the mixing angles are |141 117j 

sin^ 26*23 > 0.92 

0.70 < sin^ 26*12 < 0.95 



sin^ 6*13 < 0.048. 



(34) 
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and nothing is known so far about the CP violating phases. The existence of the two large 
mixing angles 612 and 823 is the most striking difference between the MNS and the CKM 
matrix. 

Both the mass matrices for charged leptons and for neutrino have to be diagonalized 
and their lefthanded unitary rotation matrices give rise to the MNS matrix. For the charged 
leptons the diagonalization procedure is analogous to the ones for the quarks: 

/ me \ 

























rrir 



m^'> = I \ =Vl''^M^''V^' (35) 



After the seesaw mechanism the neutrino mass matrix is 

Mf™) = M(^)-i^M(^)^ (36) 

where all ingredients are 3x3 matrices. Note that Mq*^'^"'^ has to be symmetric. The Dirac 
part is diagonalized as follows: 

1 

m(^)=| 0' mo,2 I = )t^(^)yf ) (37) 
mD,3 

Plugging the diagonalization of the Dirac part into Eq. I|36|) yields 
where the central matrix C is defined as 



Mq 

The central matrix C is diagonalized with a unitary matrix 

j mi Q \ 
C = T my = T\ m2 (40) 



mi 











m2 














The masses in the diagonal matrix m^, are the physical neutrino masses. 

The MNS matrix involves the rotations that diagonalize the mass matrices of the charged 
leptons and the neutrinos. This also includes the rotation that diagonalizes the central 
matrix. Therefore, in terms of the quantities defined above, the MNS matrix becomes 

Vmns = ^i'^V^^V. (41) 

One might hope that the striking features of the neutrino mixing matrix might allow 
us to gain some insight into the Major ana sector of the theory, which comes from very 
high scale physics and for which we cannot argue for any preferred weight. Unfortunately 
this does not happen. We have simulated the neutrino mixing by considering a variety of 
possibilities for the Majorana mass matrix. We have allowed for different power law weights 
for the Majorana mass matrix as well as flat distributions or even a common mass. In each 
case we took the Dirac Yukawa couplings to be distributed with the weight ~ 1/m^'^^ . The 
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results for the mixing angles in each case were so similar that we do not bother to display 
the differences. In each case, the mixing angles were dominated by small mixings, as we 
also found for the CKM mixing angles. 

The reason for that is that the only difference within our scheme between the CKM and 
the MNS matrix is the additional factor of the rotation matrix T that contributes to the 
MNS matrix. T diagonalizes the central matrix C from Eq. H39|) that contains two factors of 
m^^\ Since the diagonal Dirac mass matrices m^^'> come from a distribution p ~ we 
usually expect the entries of m^^^ to have a hierarchical pattern with ttid,! *C 'm£)^2 ^ itT'D,3- 
In order to obtain large mixing angles from some of the entries in the central matrix C must 
be of the same order. That however can only be achieved through a numerical coincidence 
or a correlated hierarchy of Mq^^""'^ designed to counterbalance the hierarchy present in 
m^^^ |16j . Our random matrices Mq^""'^ are always simulated independently from all Dirac 
mass terms and thus there is no correlated hierarchy and not many more cases of two large 
mixing angles than for the CKM matrix, no matter how we simulate the Majorana mass 
matrix. 

From now on, we always use Majorana mass matrices with elements distributed uni- 
formly between zero and a Majorana scale Anjaj- The charged lepton mass matrices and 
the Dirac mass matrices are simulated using a weight p ~ The neutrino masses 

we get from these simulations are close to a scale invariant distribution similar to the case 
of the quarks as seen in Fig. ^lin Section 4. 

We are interested in the probability of finding large mixing angles, especially the "two 
large mixing angle" solution found in nature, using the parameterization of Eq. (|33j) . We 
will define a large mixing angle to be an angle between 30° and 60° which corresponds to 
sin^ 26 > 0.75. 

In studying the properties of neutrino mixing, we have run simulations under various 
conditions. In the first setting, we impose no constraints at all and compare the percentage 
of two large mixing angles between the CKM matrix and the MNS matrix. We find two large 
mixing angles in 5% of the simulations for the CKM matrix, whereas for MNS matrix the 
percentage is a little higher being 6.5%. This difference is due to the seesaw mechanism i.e. 
due the effect of matrix in Eq. ()41() . Note that the resulting distributions of mixing angles 
for the MNS matrix are independent of the Majorana scale when there are no constraints. 

We also have run simulations with constraints on the charged lepton masses and the 
neutrino masses. For the charged lepton masses we adopt similar constraints as we did in 
the previous section for the quarks considering only those mass matrices for which the mass 
eigenvalues for the electron (i.e the lightest eigenvalue of the charged lepton matrix) are 
less than 5 MeV, the mass eigenvalues for the muon are within 25 MeV and 1 GeV and the 
mass eigenvalues for the tau are greater than 1.5 GeV. For the neutrino masses we explore 
two different kinds of conditions. In one set of constraints, we accept only those neutrino 
masses where the mass differences lie within the experimental 3a allowed ranges and the 
cosmological bound is satisfied. A second set of simulations only requires that the two 
heaviest masses lie in the two decade wide mass range 0.001 eV < m2, < 0.1 eV with a 
mass ratio 7713/7712 < 10, where 7773 and 7772 are the heaviest and second heaviest neutrinos 
respectively. This is more general than the first set of neutrino mass constraints, but 
captures what could be important features of the physical masses. In both the constrained 
cases, the number of cases with two large mixing angles depends on the Majorana scale. 

In Tables 0] and El we present the results for the simulations concerning the number of 
large mixing angles when the above constraints are imposed. We see that the probability 
of two large mixing angles is never dominant. However, this probability is not small either. 
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Table 4: Percentages of outcomes depending on the number of large mixing angles and the Majorana 
scale when the masses satisfy the 3a bounds for the neutrino mass differences and the cosmological 
bound. 
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43.2% 
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5.0% 



Table 5: Percentages of outcomes depending on the number of large mixing angles and the Majorana 
scale when the masses m2 and are in the range 0.001 — 0.1 eV with a mass ratio m^/m2 < 10. 



Depending on the scenario considered and the Majorana scale, the percentage ranges from 
about 4% to 18% as we vary Aj^aj- There is a general trend in the cases where constraints 
on the masses are imposed - the percentage is larger when the Majorana scale is smaller. 
This is readily understandable. When the Majorana scale is small, the average element in 
the mass matrix is comparable to or greater than the resulting mass eigenvalues, because 
in the seesaw mechanism the masses are proportional to 1/h.Maj- When all elements are 
large, one requires larger mixing angles to diagonalize the mass matrix. 

Our results are somewhat inconclusive. One finds solutions with two large mixing angles 
with a reasonable percentage. However, it is not the most likely outcome. Statistically, we 
can not conclude anything significant from a comparison with the observed mixing angles. 
This conclusion is independent of how we have modeled the possible randomness of the 
Majorana sector. We have also found that the information on the neutrino mass differences 
is unable to tightly constrain the scale of the Majorana sector. 

We can also extract some predictive elements from our analysis. We will study these 
using a Majorana scale of A^aj = 10^^ GeV. In each case we also only consider those 
realistic solutions that satisfy the 3cj bounds on the neutrino mass differences as well as the 
cosmological bound and that have large angles 9i2 and ^23 and values of ^13 smaller than 
the experimental bound, and thus are very similar to the real data. 

Using these realistic simulations, we have studied the size of the third mixing angle 
613. The distribution of sin^is is shown in Fig. within its 3cj allowed range, that is for 
sin^ ^13 < 0.048. We see that the distribution of sin^is is not peaked at low values but 
rises towards its 3(t upper bound, so that there is no reason to expect that the third mixing 
angle is suppressed. The median of the distribution is 0.14 and 95% of the values lie above 
0.04. 

We also have studied the neutrino equivalent of the Jarlskog invariant that signals the 
strength of CP violation for lepton number conserving processes the neutrino sector. This 
is plotted in Fig. 1161 Here we see that the Jarlskog invariant is generally large, of order 
10~^. We see that the peak of the distribution is close to its theoretical maximum possible 
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Figure 15: The distribution of sin^ia 



value which is 1/(6^/3). The median of the distribution is at 0.016 and 95% of the values 
are above 1.6 x 10~^. 

Finally we consider the effective Major ana mass defined by 



which is the figure of merit in neutrinoless double beta decay j20j . The simulations of this 
parameter are shown in Fig. El again for the realistic case very similar to nature. We see 
that the distribution decreases approximately linearly until a value of 0.008 eV, with only 
very few values higher than that. These come from cases where the neutrino masses are not 
strictly hierarchical and the lightest mass is close to the second lightest mass. The median 
of the distribution is at 0.0028 eV, 95% are higher than 3.5 x 10""^ eV. 95% of the values lie 
below 0.0068 eV which is consistent with other results that expect the effective Majorana 
mass for the case of hierarchical neutrino masses to be smaller than 0.0064 eV [20] . 

A pioneering effort in simulating neutrino mixing angles randomly in the context of the 
seesaw mechanism was the work of Goldman and Stephenson in 1981 jl8j . These authors 
randomly populated both the Dirac and Majorana mass matrices with a uniform distribution 
and found that the typical case involved angles of the same size as the Cabibbo angle. The 
intent was not to invoke a random dynamical framework but to argue that unless the initial 
Yukawa couplings were peculiar, that we should expect to find Cabibbo-like mixing angles. 
More recently, the idea has surfaced under the label of "anarchy" in the works of |19j . 
These authors build in more of what is now known about the neutrino mass spectrum and, 
although again one most typically finds Cabibbo-like angles, they point out that one finds 
two large mixing angles (defined as sin^ 26 > 0.5) in about 13% of the cases when simulated 




(42) 
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Figure 16: Distribution of the Jarlskog invariant J for neutrinos 

with a uniform random distribution. Our considerations have been somewhat different. We 
have used the Dirac component of the neutrino masses to be the same distribution that has 
been used in the description of the charged leptons and the quarks, and considered various 
possibihties for the Majorana sector. However, our numerical results are fairly similar. Two 
large mixing angles are found between 5% and 15% of the time as we change h-Maj and when 
we consider output masses similar to those needed experimentally. Note that the definition 
of "large mixing angles" is more generous in the anarchy papers (sin^ 20 > 0.5) than in 
our work {sm^2d > 0.75). In our simulations this can make a factor of two difference 
in the percentage of cases that are said to have two large mixing angles. For example, 
with AMaj = 10^^ GeV, the percentage of two large mixing angle solutions goes from 5.4% 
to 10.1% when adopting the more generous range. However, this does not change our 
conclusions. Large mixing angles are found a significant fraction of the time, but are not 
the most probable solution. 

Our work has concentrated on the standard seesaw mechanism. If considering the Type 
II seesaw or triplet seesaw mechanisms instead, the neutrino mass matrix is com- 
pletely determined by physics at a very high scale for which we have no basis to prefer a 
certain weight. This leads to more freedom because of the lack of constraints, but unfortu- 
nately also to no predictive power. 

6 Summary and discussion 

We have provided an exploration of the distribution of quark and lepton masses and mixings, 
treating these as random variables distributed with respect to some probability distribution 
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function, which we called the weight. We have shown that this weight is close to being scale 
invariant in form, i.e p{m) ~ 1/m. In particular a statistical analysis of power law weights 
p ~ tightly constrains the power to be J = 1.02 ± 0.08. Such weights commonly 

produce small mixing angles and a hierarchy of CKM elements. Our exploration of the 
neutrino sector was somewhat inconclusive. The use of the usual seesaw mechanism, com- 
bined with the expected weight for the Dirac masses, could marginally accommodate the 
observed mass differences and mixing angles, but the existence of large mixing angles was 
not a common occurrence in this framework. However, the neutrino sector does allow us 
to make statistical predictions of this framework. A normal hierarchy of neutrino masses is 
highly favored by the scale invariant weights. We saw that the leptonic Jarlskog invariant is 
typically of order 10~^ and with 95% confidence greater than 1.6x 10~^. The third presently 
unmeasured MNS mixing angle sin^is is typically of order 0.1 and with 95% confidence it 
is larger than 0.04. The effective Majorana mass ruee is typically of order 0.001 eV. 

We are aware of many caveats about this procedure, and there axe probably others that 
we do not enumerate. Some of these issues are: 

1) Differences between quarks and leptons: The information that we have about the 
masses exists at low energy. However, it is likely that the fundamental input consists of 
parameters defined at a high energy scale. We explored the effect of this using only the 
Standard Model to run the Yukawa couplings to the GUT scale, and this did not modify 
any of our conclusions. However, wc clearly do not know the physics at intermediate scales. 
If the intermediate physics only produces logarithmic running, the residual effect is also 
probably not large. 

2) Generation structure: The assumption of statistical independence of all masses may 
fail in some obvious ways. The most suspect aspect is that we do not assume any correlation 
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between members of the same generation. In some grand unified theories, this is incorrect, 
leading for example to a relation between the r and b quark masses [2^. Perhaps in the 
fundamental theory the random distribution occurs separately for the average value of a 
generation's masses and for the ratio of masses. This pattern is distinct from that which 
we have assumed. 

3) Anthropic constraints: In a landscape picture there are inevitable anthropic con- 
straints. Some parameter sets that occur in the landscape do not allow the existence of life, 
and hence we could not find ourselves in a portion of the multiverse where these parame- 
ters occurred. These constraints may distort the mass distribution. We have provided an 
estimate of this effect in the discussion of quark and lepton masses, and it appears not to 
change our qualitative conclusions. We have not considered a possible anthropic constraint 
on neutrino masses |25j . which would disfavor masses above 1 eV. In general, it is hard to 
fully explore the effect of anthropic constraints without having the full fundamental theory. 

4) Other flavor symmetries: Many of the attempts to understand quark and lepton 
masses have focussed on adding extra flavor symmetries to the Standard Model. Our basic 
assumption asserts that vacua with these extra symmetries do not constitute a large part 
of the string theory landscape of vacua. 

From these comments, it is clear that our investigation is a rather preliminary attempt 
at phenomenology in the context of the landscape. In the end, the only fully compelling 
procedure is to identify the fundamental theory, such as string theory, that leads to the 
landscape and fully solve it. With this solution in hand one can investigate directly those 
vacua that lead to the Standard Model and calculate the statistical features of those so- 
lutions. Then one can directly address the constraints and symmetries that occur due to 
intermediate scale physics and can assess whether the fundamental theory is statistically 
compatible with the data. However, such a "top down" solution is clearly very difficult and 
we are far from being able to accomplish it. In the meantime we need to tentatively explore 
the phenomenology from the "bottom up" as best we are able. 

It is hoped that eventually this weight can be calculated from a more fundamental theory 
such as string theory. An illustrative example was given in Ref. using the Intersecting 
Brane Worlds (IBW) construction of the Standard Model There, the Yukawa couplings 
arise from the a non-perturbative effect spanning the area between the intersections of three 
pairs of branes (Aij/.), with the Yukawa interactions being exponentially suppressed in that 
area, 

Tijk = roe-^ (43) 

up to a phase. In this case a flat distribution in the area p{A) ~ constant yields the scale 
invariant weight in masses. The range in the allowed values of the area translates into a 
finite range in the masses. It is possible that this exponential behavior could be a more 
general feature of the landscape. Of course, this relation does not fully explain the nearly 
scale invariant distribution in the masses, but only transfers the problem to the string scale 
to understand whether the distribution in the area is nearly flat. It would be interesting to 
explore possible dynamics that could produce an ensemble of Standard Model string vacua 
and assess the resulting effect on the distribution of masses. 
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Appendix: The weight for Am^ and the Majorana scale 

When we consider the neutrino masses by themselves without mixing, the hght neutrino 
masses m,- are of the form 



m 



mi = -7^ (44) 
Mm ^ ' 

where we take Mm as a common scale. Using the scale invariant weight 

P {mD,i) = — {mD,i - mow) {m - mD,i) (45) 

'"'low ' 

for the Dirac masses, we derive the weight for the neutrino mass differences ^nif- = mj —mf 
and use the physical data for a likelihood fit of the Majorana scale Mm- 

As a first step let we calculate the weight for the neutrino masses rrii with the scale 
invariant random Dirac masses rriD %■ 



p (mi) = J dmD,i p {mD,i) 5 \^-^^ - m j 

21°g;^ V ' Mm] \Mm 



Mm 

2 

rrir, 



^ Q{mi-ML)Q{Mu-mi) (46) 



where we defined M^ = and Mu = for simplicity. We see that the neutrino 

masses are also distributed with respect to a scale invariant weight. The lower and upper 
bounds are of course changed as we would expect and the result in Eq. (|i6|) is automatically 
normalized since we started from a normalized weight. 

As a next step we calculate the weight for the neutrino masses squared, p {m 

i)- The 

calculation is analogous to the one for p {nii) in Eq. and the result is 

P i^l) = A {m^ - Ml) (m2 - m^) . (47) 

log if 

With that we will calculate the weight for Amf^ = — mf: 

p (Amfj) = / dmf / dnij p (mf) p (m^) 5 [nij — mf — Arnfj) 



[m^ - |Am2.|) (^Ml + I Am! 



log^ TTTTTs (48) 



Jlog^f lAmf^-l MfjMl 



As we would expect, pyAm~jj = p^^— Am^^j = /5^|Am^j|j, meaning that there is no 
preferred sign for Amfj. The weight for Amf- is basically proportional to one over Am?^ 
up to some logarithmic corrections. For the complete weight p ( Am?^ j , we need to include 
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the range of possible values of Am?- with a Theta function: 



hp 



2 |Am2.| 



log 



|Am|| 



Ml + \/\mf- 



-Ml-\Aml\) 



(49) 



Again, we find that our result is automatically normalized. Recall that we defined Ml 



Mm 



and Mjj = so that p ( Am: 



depends on the Majorana scale Mm as well as on m/o^ 
agrees exactly with the result of corresponding 



and m*. Our analytic result for p yAmfj 
simulations. 

One benefit of explicitly knowing the weight of an observable is that one can perform a 
likelihood fit with it. If we used the weight for Am'^ along with the physical data for Amfg 
and Amfg to fit the Majorana scale Mm using a likelihood function L = p (Am^g) P (^"^13) 1 
we would make a mistake by assuming the two measured mass differences are independent of 
each other. Therefore we have to first calculate the weight for two neutrino mass differences 
p ( Am?^-, Am?^ I that takes the correlation into account: 



p(Am,2.,Am^fc) 



dmf / drrij / dmi p (m^) p (m|) p (m|) 



M^ 



1 



■ / dm^ — ^ 

J mj 



5 (m| - mf - Am^_,) 5 {ml - - Amf^ 



mf + Am^. mf + Am^^ 



e(m2-M|)G(M2 
e K - (Ml - Ami)) ii^u - Amf,) - m^) 



e m- 



1 



Am?- Am 
'J 



log 



{Ml 



\ml)) G {{M, 



mj) 



Am^ 



+ 



Am2. (Amf^. 



log 



+ 



Ml + Am 



Am^^, [Amff^ - Am 

Q{Aml)Q{{Ml-Ml)-Aml) 
e{Aml)e{Aml-Am^j) + ... 



log 

^ Ml + Ami 



(50) 



Here we only display the part of the result for the case where < Am?- < Am% which is 



ik 

the case for a normal neutrino mass hierarchy when we identify Am?- with Am^2 and Am,^ 



ik 



with Am^g. Again, this weight depends on the Majorana scale Mm via M^ and Mjj. 

Now we can finally perform a likelihood fit for the Majorana scale Mm using the like- 
lihood function L {Mm) = p (Am^gi Am^g). As the experimental input we use the cental 
values 112] Amf2 = 8.0 x 10"^ eV^ and Amfg = 2.3 x IQ-^ eV^. We find that Mm has to 
be smaller than 10^^ GeV to accommodate the experimental data. As a preferred value the 
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fit yields 9 x 10^^ GeV. At the la level Mm < 3 x 10^° GeV are excluded and at the 2c7 
level Mm < 10^ GeV are excluded. 

To interpret these results it is useful to have a second look at Fig. where we assume 
hierarchical neutrino masses mi <C m2 <C m^. For high Mm there is more "phase space" 
for the lightest neutrino mass which explains the preferred value close to its upper limit. 
But since this "phase space" decreases only linearly when one goes to smaller Mm on 
a logarithmic scale the the range of allowed Majorana scales Mm spans many orders of 
magnitude. 
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